Thermal Casimir Effect in the Plane-Sphere Geometry 
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The thermal Casimir force between two metallic plates is known to depend on the description 
of material properties. For large separations the dissipative Drude model leads to a force a factor 
of 2 smaller than the lossless plasma model. Here we show that the plane-sphere geometry, in 
which current experiment are performed, decreases this ratio to a factor of 3/2, as revealed by 
exact numerical and large distance analytical calculations. For perfect reflectors, we find a repulsive 
contribution of thermal photons to the force and negative entropy values at intermediate distances. 



Interest in the Casimir force, arising due to the con- 
finement of the vacuum fluctuations of electromagnetic 
fields between two reflecting bodies, has been continu- 
ously growing during the last decade and the motiva- 
tion for measuring it precisely has led to a number of 
original experiments using various modern technologies 
[Ell SHUSH- The Casimir force depends on a num- 
ber of factors, including the bodies' material properties 
[1, [§], their surface state [HI EE E3 an d shape. Cur- 
rent experiments are performed using a spherical and a 
plane surface. The force in this plane-sphere geometry 
is usually calculated within the Proximity Force Approx- 
imation (PFA) [l3[ which averages the force calculated 
in the plane-plane geometry over the local inter-surface 
distances. Only recently studies have been devoted to 
Casimir force evaluations goin g beyond the domain of 
validity of PFA 0, El EE EE El, El • 

The influence of temperature on the Casimir effect 
has given rise to intense discussions over the last decade 
[20l l2l| , in particular because the force exhibits an un- 
expectedly strong correlation with the detailed descrip- 
tion of optical properties of the metallic surfaces used 
in the experiments. The dielectric function of metals is 
often modeled by the plasma model where the plasma 
frequency wp, depending on the density of conduction 
electrons, acts as a cut-off frequency above which the re- 
flectivity goes to zero. In the Drudc model which also 
accounts for the relaxation of conduction electrons, the 

dielectric function at imaginary frequencies u> = it; is 

2 

given by = 1 + g7|j?^ ; with the relaxation frequency 

2 

7 related to the reduced dc conductivity a — ^ (see 
[11). The plasma model permittivity is obtained from 
the Drude model one in the limit 7^0, while the per- 
fect reflector's infinite permittivity is recovered with the 
further limit top — > 00. 

While the dissipative Drude model seems a more ap- 
propriate description of metallic mirrors, it turns out that 
recent force measurements are in better agreement with 
the predictions of the lossless plasma model [l| . In addi- 
tion, the Casimir force between two plates at large sep- 
arations turns out to be a factor of two smaller when 
calculated with the Drude model compared to the one 



obtained with the plasma or the perfect reflector model. 
This significant difference is attributed to the vanishing 
contribution of TE modes at zero frequency for dissipa- 
tive mirrors entailing that for the Casimir force, contrary 
to the dielectric function, there is no continuity from the 
Drude to the plasma model at the limit of a vanishing 
relaxation [221 ]. In contrast to the other two models, the 
Drudc model also leads to negative Casimir entropy val- 
ues between two plates [231 ] . 

In the present letter, we treat plane and spherical 
metallic surfaces coupled to electromagnetic vacuum and 
thermal fields with material properties described by ei- 
ther the perfect reflector, plasma or Drude models, and 
show that the above mentioned features are considerably 
altered in this situation. First, the factor of 2 between 
the long distance forces in Drudc and plasma models is 
reduced to a factor of 3/2 decreasing even more below 
this value when small spheres are considered. Second, 
negative entropies are found also for the perfect reflec- 
tor model, in which case they can only be related to the 
plane-sphere geometry and not to dissipation. Finally, 
PFA underestimates the Casimir force within the Drudc 
model for short distances, while it overestimates it at all 
distances for the perfect reflector and plasma model. 

We consider a large parameter space generated by the 
five length scales involved in the problem : the surface 
separation L, the sphere radius R, the thermal wave- 
length At = hc/ksT at temperature T, the plasma wave- 
length Ap = 2itc/ujy> and the wavelength associated with 
relaxation frequency A 7 = 2nc/"f. The general scatter- 
ing formula [24j for the Casimir free energy T between a 
plane and a sphere is given by: 

^ = fc B T£logdet(l-A4(£ n )) An = 2 -^^- 

n 

M(£ n ) = ^ s (e„)e-' C( «" )£ ^p(^)e-' c ^) £ . (1) 

According to ([1]), the operator M. contains the reflec- 
tion operators 7?-s aud 7?.p of the sphere and the plane 
respectively, the latter being evaluated with reference 
points placed at the sphere center and at its projection on 
the plane, as well as the translation operators e~ K ^ n%>c 
describing one-way propagation between the reference 
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FIG. 1: Thermal Casimir force at T = 300K computed be- 
tween perfectly reflecting sphere and plane, divided by the 
zero temperature force. Solid lines from bottom to top corre- 
spond to increasing values of sphere radii. The upper dashed 
curve represents the PFA expression while the lower dotted 
curve is the analytical asymptotic expression in the L S> R 
limit [Colors online]. 



points on a distance C = L + R; the primed sum is a 
sum over the Matsubara frequencies (n > 0) with the 
n = term counted for a half. 

The upper expression is conveniently written through 
a decomposition on suitable plane-wave and multipole 
bases. The resulting elements of the matrix IZp are pro- 
portional to the Fresnel reflection coefficients r p with 
p =TE and TM for the two electromagnetic polariza- 
tions, while those of IZs are proportional to the Mie co- 
efficients ag, bg [25j for electric and magnetic multipoles at 
order £= 1,2,..., respectively. Due to rotational symme- 
try around the z-axis, each eigenvalue of the angular mo- 
mentum m gives a separate contribution to the Casimir 
free energy T 1 -" 1 ^, obtained through the same formula as 
([T]), with M. reduced to the block matrix collect- 
ing the couplings for a fixed value of m. The numerical 
computations presented below are done after truncating 
the vector space at some maximum value ^ max of the or- 
bital number I. The results are thus accurate only for 
R/L smaller than some value which increases with £ max , 
typically R/L < 10 for our £ max = 34, remaining a factor 
of 10 below current experimental values R/L > 10 2 . 

The results of the numerical computations are shown 
on Fig.l for perfect reflectors. We have calculated the 
Casimir force F perl between the plane and the sphere at 
ambient temperature and then plotted the ratio $ perf of 
this force to its value at zero temperature: 



F pcrl (L, T) = - 



dL 



erf 



FP crf (L,T) 
FP CTf (i,0) ' 



(2) 



The various solid curves are drawn for different radii R of 
the sphere as a function of the distance L, with increasing 
values of R from bottom to top ; the upper dashed curve 
on Fig. 1 represents the quantity $ppA as ^ would be ob- 



FIG. 2: Same plot as in Fig.l for the Casimir force at T = 
300K computed with the Drude model, divided by its value at 
zero temperature. The dashed curve correspond to the PFA 
expression [Colors online]. 



tained from ([2]) by using PFA ; the lower dotted curve is 
an analytical asymptotic expression discussed below. 

Fig. 2 shows the variation of the ratio $ Drud ; defined as 
in ^ for the Drude model with Ap = 136nm, A 7 /Ap = 
250 and At = 7.6fim. The dashed curve on Fig. 2 rep- 
resents ^ppA d as obtained for the Drude model by using 
PFA. We do not plot the variation of ^ plas , defined as 
in ([2|) for the plasma model, since it is as expected quite 
close to the one shown on Fig.l for perfect mirrors. 

In most cases, the ratio 1?, starting from unity at small 
distances, decreases below unity when the distance in- 
creases, then reaches a minimum before increasing at 
very large distances. While such behavior was already 
observed for dissipative Drude mirrors in the plane-plane 
geometry [U [23| (the dashed PFA curve on Fig. 2 is also 
below unity for L < 4/zm) , or in the atom-surface config- 
uration out of equilibrium [26|, our computations show 
quite unexpectedly that in the sphere-plane geometry 
such a behavior takes place even for the perfect reflec- 
tor and plasma models at thermodynamical equilibrium. 
Hence, for all three models the contribution of the ther- 
mal photons to the Casimir force can be repulsive, which 
suggests that the entropy could be negative for some val- 
ues of the parameters (see below). 

A second important feature comes out in a striking 
manner from the comparison of Fig.l and Fig. 2 : the 
PFA always overestimates the effect of temperature on 
the force between mirrors described by a perfect reflector 
or plasma model; in contrast, it underestimates the tem- 
perature effect for the Drude model at small distances, 
while it overestimates it at large distances. The over- 
estimation is however smaller than for perfect mirrors. 
As a consequence, for small separations the Drude and 
plasma models lead to Casimir force values much closer 
than predicted by PFA. These results clearly demonstrate 
the strong correlation between the effects of plane-sphere 
geometry, temperature and dissipation. 

In the following we will corroborate the previous nu- 
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merical results by presenting analytical calculations of 
the thermal Casimir energy in the limit of large distances 
(L ^> R). Since the number of modes € max needed to 
get an accurate result decreases when L/R, increases, we 
may take ^ max = 1 in this limit. Another consequence of 
this limit is that the reduced frequency £ = £i?/c is very 
small, since the characteristic frequencies, which give the 
main contribution to the Casimir force, scale as £ ~ c/L. 

For perfect reflectors, where Ap = 0, the dielectric 
function e is infinite at all frequencies and we obtain 
the following low-frequency expansions for the Mie coef- 
ficients a\ and b\ describing the scattering on the sphere: 



pcrf 
i — 



(3) 



2|3 

3 ' _1 3 

The other steps in the calculation of the Casimir force 
may then be done analytically and the sum over all Mat- 
subara frequencies may be given in a closed form. One 
obtains in this manner the following approximation of 
the Casimir free energy 



4A T £ 3 ' A T 

_ v sinh v + cosh v{v 2 + sinh 2 v) 
2 sinh 3 v 



The fact that the upper expression is a relevant approx- 
imation is shown on Fig.l: the lower dotted curve, rep- 
resenting the value of the ratio $ perf deduced as in (J2|) 
through a derivation of expression (|4]), is close to curves 
computed for small radii R <Si L. Using (j4]), it is straight- 
forward to derive an analytical expression of the entropy 
S = -dTjdT: 

S perf =^^(<^)+^>>)), L»R , (5) 

which gives negative values for v < 1.5, that is L < 1.8/im 
at T = 300K. 

In addition, we can derive from ^ low- and high- 
temperature expressions for the free energy: 



erf 



erf 



9hcR 3 
' 16ttL 4 
3hcR 3 



1 



' 135 + 945 
L > A T , R- 



A t > L > R, 



8XtL3 , (6) 

Equation ^ is the large-distance high-temperature limit 
which can be generalized to metallic scatterers described 
by either the plasma or the Drude model. Starting with 
the lossless plasma model (7 = 0) we obtain for L ^ Ap 
Fresnel coefficients with unit modulus tte ~ — 1, ?"tm ~ 
1, while the low- frequency expansion of the Mie coeffi- 
cients [13], and the resulting free energy, are read, intro- 
ducing the parameter a = ^f^-'- 



plas 



2£ LPlas 



1 



coth< 



3hcR 3 ( 1 
~8A T £ 3 \ 1+ 



cotha 
a 



e 3 , (7) 



, L > At, R- 



The result for perfect reflection is reproduced by ([7]) when 
both L, R > A P . 

For the dissipative Drude model (7 ^ 0), the low fre- 
quency limit of the two Fresnel coefficients have the well- 
known form tte - * 0, ttm ~ 1- The low- frequency ex- 
pansion of the Mie coefficients and the ensuing free en- 
ergy are read 



Drud 

"1 



Drud 



2_l 

3 ' (TqR 
hcR 3 



tDrud 



45c 



4A T L 3 



L > At, R- 



(8) 



The long distance free energy for the Drude model 
amounts to 2/3 of the value for perfect mirrors whereas 
this ratio is 1/2 in the plane-plane geometry. The latter 
result is explained by the fact that the TE reflection co- 
efficient vanishes at zero frequency so that only the TM 
modes contribute [2(J[2l|]. The change of the ratio 1/2 to 
2/3 in the plane-sphere geometry has to be attributed to 
the redistribution of the TE and TM contributions into 
electric and magnetic spherical eigenmodes. 

Formally the results for the Drude model © can be 
obtained from the plasma model results by taking the 
limit R <C Ap. In this limit however, we should take into 
account the effect of quantum confinement in the small 
sphere, which is out of the scope of the present letter. 
Two further features in ([5]) must be emphasized. First, 
the coefficient b\ is vanishingly small in the Drude model 
but not in the plasma model; the latter can thus not be 
obtained by turning the relaxation frequency 7 to zero 
(or (To to 00). In addition, the free energy for the Drude 
model is independent of the values of Ap and 7, whereas 
the one for the plasma model depends on Ap. 

On Fig. 3, we illustrate the comparison of the two mod- 
els by plotting the ratio of the thermal Casimir forces 
F plas calculated with the plasma model and F Dlud ob- 
tained with the Drude model. Again, the plots corre- 
spond to Ap = 136nm and A 7 /Ap = 250. The results 
of our calculations are shown by the solid curves with 
the sphere radius increasing from bottom to top as in 
Fig.l. The ratio jTP las /^ Drud varies in the plane-sphere 
geometry as a function of the sphere radius, which clearly 
demonstrates the strong interplay between the effects of 
temperature, dissipation and geometry. For large spheres 
(R S> Ap), the ratio converges at long distances to the 
value 3/2 which has been obtained analytically in the pre- 
ceding paragraphs, whereas it remains smaller for small 
spheres (down to 1.2 for R ~ lOOnm). The dashed curve 
gives the variation of the same ratio as calculated within 
the PFA which leads to a factor of 2 in the limits of large 
distances or high temperatures. We emphasize that the 
factor of 2 deduced within PFA is never approached at 
the large distance limit within the calculations performed 
in the plane-sphere geometry. 

To summarize we have computed exact results for the 
Casimir free energy and force at non zero temperature 
in the plane-sphere geometry. We have given plain evi- 
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in the relative approaching of the Casimir force values 
computed with the Drude and plasma model, the ap- 
pearance of negative entropies evidently not related to 
the presence of dissipation and the fact that PFA under- 
estimates the Casimir force for the Drude model at short 
distances while it overestimates it for the plasma model. 
If the latter feature were conserved for the experimental 
parameter region R/L (> 10 2 ), the actual values of the 
Casimir force calculated within plasma and Drude model 
could turn out to be closer than what PFA suggests, even- 
tually diminishing the discrepancy between experimental 
results and predictions of the thermal Casimir force using 
the Drude model. Settling this question is an open and 
highly topical program in Casimir physics. 



FIG. 3: Ratio of thermal Casimir force at T = 300K calcu- 
lated with the plasma model and the Drude model, as a func- 
tion of surface separation L for different radii of the sphere. 
The solid curves from bottom to top correspond to increasing 
values of sphere radii. The dashed curve is the PFA predic- 
tion. [Colors online]. 

dence for a strong correlation between the effects of ge- 
ometry, temperature and dissipation based on the per- 
fect reflector, plasma and Drude model to describe ma- 
terial properties. The correlation becomes clearly visible 
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